Two-Photon Spiral Imaging with Correlated Orbital Angular Momentum States 
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The concept of correlated two-photon spiral imaging is introduced. We begin by analyzing the joint 
orbital angular momentum (OAM) spectrum of correlated photon pairs. The mutual information 
carried by the photon pairs is evaluated, and it is shown that when an object is placed in one of the 
beam paths the value of the mutual information is strongly dependent on object shape and is closely 
related to the degree of rotational symmetry present. After analyzing the effect of the object on the 
OAM correlations, the method of correlated spiral imaging is described. We first present a version 
using parametric downconversion, in which entangled pairs of photons with opposite OAM values 
are produced, placing an object in the path of one beam. We then present a classical (correlated, 
but non-entangled) version. The relative problems and benefits of the classical versus entangled 
configurations are discussed. The prospect is raised of carrying out compressive imaging via two- 
photon OAM detection to reconstruct sparse objects with few measurements. 

PACS numbers: 42.30.Va,42.50.Tx,42.65.Lm 
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I. INTRODUCTION 

In digital spiral imaging (DSI) [H, an object is il- 
luminated by light with a known spatial distribution, 
or equivalently, of known orbital angular momentum 
(OAM) distribution, and the OAM spectrum after 
the object is measured. The shape of the outgoing spec- 
trum allows determination of some properties of the ob- 
ject, or possibly identification of the object from a known 
set. However, as we will see below, this method is inca- 
pable of reconstructing the actual shape of the object; de- 
spite its name, it is inherently a non-imaging technique. 

Here, we propose correlated spiral imaging (CSI), mea- 
suring correlations of OAM values within two-photon 
states or between two light beams. We will consider 
measurement of the correlations through both coinci- 
dence counting and through interference between the two 
beams. We will then show that (i) the CSI coincidence 
rate displays clear signatures of object spatial properties, 
(ii) the mutual information carried by the detected pair 
has a strong dependence on object shape and measures 
the object's rotational symmetry, and (iii) a version of the 
setup does allow efficient reconstruction of object shape, 
opening up the possibility of carrying out compressive 
imaging with high-dimensional OAM states. 

The experiments proposed here differ significantly from 
that carried out in In the latter, after filtering for 
specific OAM values, the spatial locations of the outgoing 
photons are measured in one arm, as in traditional ghost 
imaging [6-9]. But in the present case, no information 
about the spatial location or momentum of the photon is 
recorded; only angular momentum values are detected. 

In the following sections, we describe two categories of 
correlated spiral imaging experiments, one involving en- 
tangled photon pairs produced via downconversion, the 
other using classically correlated beams. The key point 



in all of the variations we describe is that there are two 
light beams (or two photons) with correlated OAM val- 
ues. Our purposes here are two-fold: both scientific and 
applied. On the pure science side, the entangled version 
is of great interest, following in a direct line from work 
such as that of [lOj and |5j, and offering a new window 
into both the downconversion process and the quantum 
correlations between the signal and idler OAM values. 
The latter correlations are certainly of scientific interest 
in their own right, apart from any applications. On the 
applied side, the classical version is likely to be more use- 
ful, as discussed more fully in section IVl 



II. BACKGROUND 

A. Laguerre-Gauss modes 

We decompose ingoing and outgoing beams in terms 
of optical Laguerre-Gauss (LG) modes. The LG wave- 
function with OAM lh and with p radial nodes is 
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and the arctangent term is the Gouy phase. 
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FIG. 1: Digital spiral imaging: the presence of an object in 
the light beam alters the distribution of angular momentum 
values in the outgoing light. 



B. Digital spiral imaging 

DSI [l| is a form of angular momentum spectroscopy 
in which properties of an object are reconstructed based 
on how it alters the OAM spectrum of light used to il- 
luminate it (fig. [1]). The input and output light may 
be expanded in LG functions, with the object acting 
by transforming the coefficients of the ingoing expansion 
into those of the outgoing expansion. Information about 
the transmission profiles of both phase and amplitude 
objects may be retrieved pi [l2|. 

The idea naturally arises of trying to use the mea- 
sured OAM spectrum to reconstruct an image of the ob- 
ject. But, although a great deal of information may be 
obtained about the object in this manner, it is not suf- 
ficient to reconstruct a full image of the transmission or 
reflection profile. To see this, expand the output am- 
plitude according to Yli p A\ v u\ v . Projecting out partic- 
ular / and p values, the detector tells us the intensity of 
each component, allowing the \ Ai p \ 2 to be found, with no 
phase information retrieved. We thus have an incoherent 
imaging setup, with total detected intensity of the form 

i p \Aip\ 2 \ui p \ 2 . But the quantities \ui p \ 2 are rotationally 
symmetric for all values of / and p (see the right-most 
panel of fig. [2]). Any image built from them is also sym- 
metric; variation of the object about the axis is lost. In 
contrast, the real and imaginary parts are not rotation- 
ally invariant (left two panels of fig. [2]), so a coherent sum 
of the form | ^Z ip Ai p ui p \ 2 allows azimuthal structure to 
be reconstructed from the interference terms. For image 
reconstruction, we thus need to obtain a coherent super- 
position of amplitudes. This can be seen in fig. [3j an 
opaque square is placed in the beam and the two expan- 
sions (coherent and incoherent) are computed, assuming 
that only the p = components are measured and keep- 
ing terms up to \l ma x\ = 15. In the left panel, where no 
phase information is assumed, the reconstructed image is 
rotationally invariant and there is no way to distinguish 
what the actual shape of the object was. In contrast, 
the coherent expansion on the right side of the figure 
produces a recognizably square output. The phase infor- 



FIG. 2: The real and imaginary parts of the Laguerre- Gauss 
function are not rotationally-invariant, in contrast to its ab- 
solute square. This is illustrated for the case of I = 1, p = 0, 
but is true generally. 



No Phase Information With Phase Information 



i 



View from side 





View from top 



FIG. 3: Incoherent (left) and coherent (right) expansions in 
Laguerre- Gauss functions, an opaque square object. In the 
former case, all variation of the object with angle around the 
axis is lost, (p max — and i>max = 15 assumed.) 



motion is vital in reconstructing the actual image shape. 



C. Entangled OAM beams 



Consider a pump beam of spatial profile $(r) = 
ui oPo (r) encountering a \ 2 nonlinear crystal, producing 
two outgoing beams via spontaneous parametric down- 
conversion (SPDC). Assume a thin crystal located at the 
beam waist (z = 0). The output is an entangled state 
[To[, with a superposition of terms of form p ' p ' , 
angular momentum conservation requiring Iq = l[ + 1' 2 . 
We will take the pump to have Iq = 0, so that the 
OAM values just after the crystal are equal and oppo- 
site: l[ = —V 2 = I. The p f 1: p 2 values are unconstrained, 
although the amplitudes drop rapidly with increasing p' 
values (see eq. (|4|) below). The output of the crystal 
may be expanded as a superposition of signal and idler 
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LG states: 



l*>= E E C l ^l\l[,p' 1 ;l' 2 ,p' 2 )5(l -l' 1 -l' 2 ), 

l'l ^2~~ 00 Pi ^2~^ 

where the coupling coefficients are given by 

C%1= |d 2 r$(r)K p ,(r)^(r)]*. (3) 

For the case of pump beam with 1$ = po = this gives 
the coefficients [12|, [l3j : 
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VgijggKj + Pi)K* + (j + m + n)! 
(pi — m)!(p2 — + m)!(Z + n)\ ml n\ 



III. JOINT OAM SPECTRA 

We now investigate the use of two beams, rather than 
one, in combination with spiral imaging. The full benefits 
of doing this will emerge in section IIV A[ In the current 
section, we focus on examination of the OAM correla- 
tions. We begin with an entangled version, where the 
light source is parametric downconversion in a nonlinear 
crystal such as /3-barium borate (BBO). Imagine an ob- 
ject in the signal beam (fig. [4]). Since OAM conservation 
holds exactly only in the paraxial case, we assume the sig- 
nal and idler are produced in collinear downconversion, 
then directed into separate branches by a beam splitter. 
(Throughout this paper we assume all beam splitters are 
50-50.) Assume perfect detectors for simplicity (imper- 
fect detectors can be accounted for by the method in 

0). 

Let P(li,pi;l2,P2) be the joint probability for detect- 
ing signal with quantum numbers h,pi and idler with 
values lii Pi- The marginal probabilities at the two detec- 
tors (probabilities for detection of a single photon, rather 
than for coincidence detection) are 



(5) 



h,P2 



Pi(h,P2) = ^2 P(h,Pi;k,P2). (6) 

h,Pi 

Then the mutual information for the pair is 

Imax Pmax 

= E E P(h,Pi;h,P2) (7) 

h,l2=lmin Pl,P2=0 

P(h,Pl]l2,P2) 



x log 2 



Ps(ll,pl)Pi(h,P2) 



The most common experimental cases are when (i) the 
values of p\ and p2 are not measured (so all possible 
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FIG. 4: Setup for analyzing object via orbital angular momen- 
tum of entangled photon pairs. 




FIG. 5: An opaque strip of width d placed in the signal path. 
The widths are (a) d — .lwo, (b) d — .9wo, (c) d — 2.5ti>o- 
The outgoing joint angular momentum spectra are plotted. As 
the width increases, the peak in the spectrum broadens, then 
(at d = wo) splits into two peaks. 



values of pi and p2 must be summed, p m ax = oo ), or 
(ii) only the p\ = p2 = modes are detected (pmax — 
0). Except when stated otherwise, we will use l m ax — 

~lmin ~ 10 and Pmax 0. 

If the transmission function for the object is T{x), the 
coincidence probabilities P(Zi,pi; I21P2) = lAL^iLl 2 have 



amplitudes 
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u VxV > x (x\,z) [u hpi (x!,z)]* T(x 1 )d 2 x 1 (9) 



where Co is a normalization constant. Here it is assumed 
that the total distance in each branch is 2z (see fig. [4j). 

That the object's size and shape affect the coincidence 
rate is easy to see. For example, fig. [5] shows the cal- 
culated spectrum when a single opaque strip of width d 
is placed in the beam. Fig. [6] shows the corresponding 
mutual information, assuming that only the p\ = P2 = 
component is detected. In both figures, we see clear ef- 
fects of changing an object parameter (the strip width). 

The central peak of the spectrum (fig. [5]) broadens as 
d increases from zero, reducing the correlation between 
li and I2] the mutual information between them thus de- 
clines, as seen in the d/wo < 1 portion of fig. [6j But 
at d/wo w 1, the central peak in /1} space bifurcates 
into two narrower peaks (right side of fig. [5j); the in- 
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FIG. 6: Mutual information versus width of opaque strip. The 
horizontal axis is in units of wo . The minimum information 
occurs at d — wq. 



formation thus goes back up as the peaks separate, as 
indeed is the case in the d/wo > 1 region of fig. [6] If we 
continue to wider the two peaks once again broaden 
and the mutual information decays gradually to zero. In 
addition, the total intensity getting past the opaque strip 
will continue to drop, so coincidence counts decay rapidly. 



IV. MUTUAL INFORMATION AND 
SYMMETRY 

Fig. [71 shows the computed mutual information for 
several simple shapes. It can be seen that / depends 
strongly on the size and shape of the object, so that for 
object identification from among a small set a compar- 
ison of the / values rather than of the full probability 
distribution may suffice. 

If the object has rotational symmetry about the pump 
axis, then its transmission function T(r) depends only on 
radial distance r, not on azimuthal angle (j). The angu- 
lar integral in eq. (|9]) is then e'^^d^ = 2tt£ M /. 
So the joint probabilities reduce to Pih^h) = f(h)^h,i 2 
(assuming pi — p2 =0). The marginal probabilities 
for each arm reduce to Pi(h) = f(h) and P2G2) = 
f(h)- The mutual information 7(^1,^2) = Si (la) where 
Si(Li) = — ^2 Li f(h) ln/(/i) is the Shannon information 
of the object arm OAM spectrum. Thus in the case of 
rotational symmetry, the second arm becomes irrelevant 
from an information standpoint. In this sense, the quan- 
tity L2) = \I(Li, L2) — S\{L\)\ is an order parame- 
ter, capable of detecting breaking of rotational symmetry. 

More generally, suppose that the object has a rota- 
tional symmetry group of order TV; i.e., it is invariant 
under (j) — > (j) + From eqs. (pQ) and ((9]) it follows that 

the coefficients must then satisfy a 1 , 11 = e^^^^aW^ , 

which implies &^ Vl = except when is integer. 

When N goes up (enlarged symmetry group), the number 
of nonzero Q 1 ^!^ goes down; with the probability concen- 
trated in a smaller number of configurations, correlations 
increase and mutual information goes up. This may be 
seen in the three right-most objects of fig. [7J for example. 



1=2.3589 1=2.7174 1=2.5704 1=2.7589 1=3.0112 

FIG. 7: The mutual information depends strongly on size and 
shape of the object. Here, the two objects on the left have 
widths 1.5wo and .2wo; all other widths are Awo. 

A. Imaging. 

The inability of DSI to produce images due to loss of 
phase information has been pointed out. Here we show 
that a variation on the entangled CSI setup can be used 
to find the expansion coefficients including phase. 

First, we note from eq. (HI) that the factors C l y l J are 

p 1 p 2 

real and positive, so the phases of the amplitudes A^£ 2 
are entirely determined by the phases of the a~}^ 1 {z). 
Note further from eqs. (pQ) and (j9j) that the only p- 
dependence in the phase of a"/ 2 ^ 1 ^) is in the factor 
e -2ipip(z) ^ w h ere i/j(z) = arctan-^. If z is much greater 
than zr, then -0 « ^ is roughly constant, so that 
e -2ipip(z) ^ (— Since we assume the outgoing p val- 
ues are pi = P2 =0, the relevant detection amplitude 
is 

^00 ~ p[0 a p[0 ^ °00 a 00 ^°10 a 10 ' 

Pi 

(10) 

due to rapid decay of the Cy q' 2 with increasing \p[ — 

Pi\. So if we break cl~, 1 q 11 into amplitude and phase, 

a ^n' Zl = r~/ 2 '^e^ p/ i ZlZ2 , then the phase is independent 
Pi u P\ 

of p[ , except for a relative minus sign between even and 
odd p[ terms, so that 

4 1 o' a =PJi,J 2 e i * ,1,2 > (11) 

where 0/ 1 / 2 is the value of ^p^hh f° r even V\-> an d pi x i 2 = 

(Coo 2M ~ c io 2M ) r o h,h is real and positive. Thus, the 

phase of Aqq 2 is the same as the phase of a p ^o for even 
p[ and differs from that of by a factor of 7r for p[ 
odd. Finding the phases of the coincidence detection am- 
plitudes A 1 qq 2 therefore suffices to determine the phases 
of all of the a^o coefficients. 

The measurement of these phases is accomplished by 
inserting a beam splitter to mix the signal and idler 
beams before detection, as in fig. [8j erasing information 
about which photon followed which path. We then count 
singles rates in the two detection stages, rather than the 
coincidence rate. If value / is detected at a given detector 
it could have arrive by two different paths, so interference 
occurs between these two possibilities. The detection am- 
plitudes in the two sets of detectors D+ and D_ involve 
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Detectors D. 




Detectors D + 



FIG. 8: A configuration allowing image reconstruction via 
phase- sensitive measurement of entangled OAM content. 

factors A + ~ (l + iafo"* 2 '* 1 ) and A_ ~ (i + 4° ~ Wl ) , 

with detection rates #± - l + |ao° "* 2 ' Zl | 2 ±2i /m a^ 2 '* 1 . 
From these counting rates, both the amplitudes and the 
relative phases of all coefficients can be found, allowing 
full image reconstruction. 

V. CLASSICAL CSI 

In recent years, it has been shown that ghost imaging 
and other "quantum" two-photon effects may be carried 
out using classically-correlated sources |l3-[l9j . It is ap- 
parent that the same is true in the case of correlated 
spiral imaging: classical OAM correlation, rather than 
entanglement, is sufficient. The essential point in the 
present case is having two spatially separated beams such 
that if the OAM detected in one beam is known, then the 
OAM reaching the object can be predicted. So all that is 
needed is strong classical correlation or anti-correlation 
between the OAM in the two arms. 

The classical analog of apparatus of fig. [8] is shown in 
fig. [9] At the left, the system is illuminated with light 
that has a broad range of OAM values (a broad spiral 
spectrum). The beam is split, with one copy passing 
through the object, and the other entering the reference 
branch. The two beams are mixed at the beam splitter, 
then the OAM content at the two detectors is measured. 
The coefficients C l p \^ 2 will no longer be given by eq. HJ 
but instead will have values determined by the proper- 
ties of the specific input beam being used. The mutual 
information between the classical beams may be defined 
just as in eq. [71 

The classical configuration of CSI has a number of 
practical advantages over the entangled version: align- 
ment issues are greatly reduced, single photon detectors 
are not needed, and much higher brightness and counting 
rate may be obtained. There is one problem that arises, 
however, which is not present in the entangled case: if a 
broad spiral spectrum is used for the illumination, then 
there is no intrinsic correlation between the OAM value 
I2 in the reference branch and the value l[ that occurs 
between the source and object. Without this correlation, 
the value of l[ is unknown and so the change in I pro- 



Detectors D. 




Broad spiral- 
spectrum input Detectors D + 
beam 



FIG. 9: A classical version of correlated spiral imaging. An 
input beam with a broad range of OAM values is split at a 
beam splitter, sending a portion through the reference branch, 
and the rest to the object. 



duced by the object is also unknown. On the other hand, 
instead of a broad spiral spectrum, we may send in sin- 
gle OAM values, one at a time, building up the OAM 
correlation function one value of l[ at a time. But this 
slows the process of image reconstruction considerably: a 
range of OAM values needs to be scanned over, one after 
another, changing a spiral phase plate or some other type 
of OAM filter multiple times in each run. In a kind of 
quantum parallelism, the entangled version can send in 
a broad range of values simultaneously, and the entan- 
gled nature of the source will automatically ensure that 
the pairs detected are of opposite initial OAM if a short 
enough coincidence time window is used. In any case, 
whether the classical or entangled version is used, two 
correlated beams are necessary in order to reconstruct 
the relative phases of the various OAM amplitudes. 



VI. CONCLUSION: POSSIBLE EXTENSIONS 
AND VARIATIONS 

Compressive imaging. Recent years have shown an 
explosion of interest in compressive sensing [20|, Hl| , in- 
cluding compressive ghost imaging [22|. The basic idea 
is that most images are very sparse when expanded in an 
appropriate basis, with the vast majority of expansion 
coefficients being very small. So if a sampling procedure 
is used that only measures the relatively small number 
of large expansion coefficients and neglects the rest, the 
image may be reconstructed from a very small number of 
measurements, often much smaller than naively expected 
from the Shannon-Nyquist theorem. 

The joint OAM spectra (such as those shown in fig. [5j) 
have been calculated for a variety of other opaque objects 
of various shapes, and in all of them it has been found 
that only a small number of the coefficients have sig- 
nificant amplitude; LG functions can therefore serve as a 
sparse basis for these shapes. It is likely that this will also 
be true of at least some classes of more complex objects. 
The possibility thus opens of compressive imaging with 
OAM states by the CSI method. Only two additional in- 
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gredients are needed: (i) Instead of taking / = 0, as we 
did in section HIH we should illuminate the crystal with a 
broad range of Zo, providing a large number of randomly 
occurring input states; this provides the large number 
of randomly chosen sampling bases needed for compres- 
sive imaging, (ii) The basis used for sensing (LG basis, 
| /,£>)) and that used to reconstruct the ima ge (po sition 
basis \r)) should have low mutual coherence [20],|2l|; this 
means the maximum value of |(r|/p)| = \ui p (r)\ should 
be as small as possible, implying that the range of I and 
p values used should be centered at the largest possible 
mean value. A detailed discussion of this will be carried 
out elsewhere. 



OAM Ghost imaging. We have seen that the pres- 
ence of an object yields clear signatures in both the CSI 
coincidence rate and the mutual information. A varia- 
tion to be explored elsewhere is "ghost" spiral imaging, 
where the l\ values are not measured (analogous to hav- 
ing a bucket detector in ordinary ghost imaging). That 
signatures of the object will still appear can be seen by 
imagining summing over the l\ rows in the histograms of 
fig- El 
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